Universality in nuclear dependence coefficient a(qr) 
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We derive the nuclear dependence coefficient a(qx) for 
Drell-Yan and J /ij) production. We show that at small qx, the 
a(qx) is given by an universal functional form: a(qx) — a + 
b q\ , and the parameters a and b are completely determined by 
either perturbatively calculable or independently measurable 
quantities. This universal functional form ct(qx) is insensitive 
to the A, and is consistent with existing data. 
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Anomalous nuclear dependence in transverse momen- 
tum spectrum of single hadron production in high energy 
hadron- nucleus collisions, known as Cronin effect jjj , was 
first observed nearly 25 years ago. In recent years, new 
data on nuclear dependence for the Drell-Yan [||^) and 
J/tp production ^-0] have renewed our interests in the 
observed novel effect However, due to technical dif- 
ficulties in handling multi-scale calculations, QCD pre- 
dictions for such nuclear dependence only exist for large 
transverse momentum (qr) region p-pd|, while almost 
all existing data are in small qx region j g-fl . In this let- 
ter, we derive the nuclear dependence coefficient a(qx) 
for Drell-Yan and J/ip production in small qx region. 

For the Drell-Yan production, the a(qx) is defined as 
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where Q is the Drell-Yan pair's total invariant mass, A 
is the atomic weight of the target, and da hA /dQ 2 dq T 
and da hN /dQ 2 dq T are the transverse momentum spec- 
trum in hadron-nucleus and hadron-nucleon collisions, 
respectively. We argue below that in small qx region 
(<7t < 2.5 GeV), the Drell-Yan g^-spectrum at fixed tar- 
get energies in both hadron-hadron and hadron-nucleus 
collisions can be well represented by a Gaussian-like dis- 
tribution, and corresponding parameters can be com- 
pletely determined by the quantities that are either cal- 
culable in QCD perturbation theory or independently 
measurable in experiments. From such a Gaussian-like 
spectrum, we derive our universal functional form for the 
a(q T ). 

Depending on physics origins, the Drell-Yan qx- 
spectrum can be divided into three regions: (I) small 
qT region where qx < qx ~ 1 GeV; (II) intermedi- 
ate qx region where q§, < qx < qx with qtj, = kQ 



and k ~ 1/3—1/2; and (III) large qx region where 
qx > qx- in region (I), the Drell-Yan g^-spectrum in 
hadron-nucleon collisions is dominated by the intrinsic 
transverse momenta of colliding partons. In perturba- 
tive calculations, the spectrum at the leading order is 
proportional to S 2 (qx)- If the partons' intrinsic trans- 
verse momentum has a Gaussian-like distribution, its ef- 
fects can be included by replacing the (5-function with 
5{q x ) = lim T _>o J- exp[— q 2 /2r 2 }, and we obtain the 



Drell-Yan gx-spectrum in region (I): 
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with Noy a dimensional normalization and r the width. 

In region (II), the spectrum can be calculated in QCD 
perturbation theory with resummation of large loga- 
rithms, such as (a s \n 2 (Q 2 / q T )) n , which are due to the 
gluon radiations from incoming partons |l3| . The resum- 
mation is extremely important for W ± and Z° produc- 
tion at collider energies because the Q 2 /q T can be as 
large as 8 x 10 3 for qx ~ 1 GeV. However, for the Drell- 
Yan production at fixed target energies, the resummation 
is much less important because of much smaller value of 
Q 2 /q T - In fact, as shown in Ref. p2] |, all existing data 
for qx as large as 2.5 GeV can be well represented by 
extending the Gaussian-like distribution in Eq. (^|) from 
region (I) to region (II). 

In region (III) , the transverse momentum spectrum can 
be calculated in perturbative QCD Q. Therefore, the 
Drell-Yan g-r-spectrum at fixed target energies can be 
represented by 
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where da^ 111 ) /dQ 2 dq\ is the perturbatively calculated 
gT-spectrum for qx > g^, and da^ /dQ 2 dq T , defined in 
Eq. (U), fits data in regions (I)+(II). 

Using moments of the Drell-Yan g^-spectrum in 
Eq. (|J), we can relate Ndy and r in Eq. (g) to physi- 
cal quantities. At fixed target energies, the contribution 
from the second term in Eq. (^) to da/dQ 2 is much less 
than one percent. Therefore, up to less than one percent 
uncertainty, Nx>y ~ da/dQ 2 . 

Define the averaged transverse momentum square as 
(g 2 ) ee / dq\ q\ (da /dQ 2 dq T )/ (da/dQ 2 ). We find || 
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that the contribution of the second term in Eq. (|^) to (gf,) 
is much less than ten percent of the first term. Therefore, 
by iteration, we obtain 2r 2 « (g^) — r(g^) with 



r(g£) 
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where 2r 2 in da^ jdQ 2 dq\ in Eq. ^ is approximately 
given by (g^). Substituting the N DY and r into Eq. (0), 
we obtain the Drell-Yan spectrum for q? < qlf as 
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Most importantly, the T(q^) in Eq. (||) is small and per- 
turbatively calculable. 

For the Drell-Yan g^-spectrum in hadron-nucleus colli- 
sions, we also need to consider multiple scattering. Sim- 
ilar to the single-scattering case, at the leading order in 
perturbation theory, the double-scattering contribution 
is also proportional to a 8- function p5[, 
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where the subscript D indicates the double scattering, 
T a „(x. , xi, X2, fcr) is the quark-gluon correlation function 
pq , [l6| , where x, X\, and X2 are the momentum frac- 
tions carried by the quark and gluon fields. The kx in 
Eq. (|^) represents the intrinsic momentum of the gluon 
which gives additional scattering. Following the same 
arguments leading to Eq. (||), we can show |l2| that if 
the partons' intrinsic /cT-dependence has a Gaussian-like 
distribution, the double scattering contributions to the 
Drell-Yan gr-spectrum in small qr region can also be 
represented by a Gaussian form. 

In high qr region, the Drell-Yan g^-spectrum in 
hadron-nucleus collision also has a perturbative tail. The 
Drell-Yan g^-spectrum at large or in hadron-nucleus col- 
lisions was calculated in Ref. flO(| . The nuclear depen- 
dence of Drell-Yan g*r-spectrum depends on two types 
of multiparton correlation functions inside the nucleus: 
T DH and T SH , which correspond to the double- hard 
and soft-hard double scattering subprocesses respectively 
p0| . These correlation functions are as fundamental 
as the well-known parton distributions, and can be ex- 
tracted from other physical observables. 

Similar to deriving Eq. (0), we can derive 



contribution to (q^) hA , and it is calculable and depends 
on the perturbative tail of the g^-spectrum. 

Substituting Eqs. (|) and @) into Eq. ([[]), we derive 
a(g-r) for the Drell-Yan production in small g-r region: 
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where AT{q^) hA = T(qk) hA -T(q^) hN , and is much 
smaller than A(g 2 )^ @. The a DY {q T ) in Eq. (§) has 
a quadratic dependence on qx- 
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At the leading order in a s , A(g^) 
For normal nuclear targets, T^ g H (x) = \ 2 A 4 ^ 3 q(x) with 
q(x) the normal quark distribution and the parameter A 2 
proportional to the size of averaged color field strength 
square inside a nuclear target || . Consequently, we have 
A(q^) hA = b DY A 1 / 3 with b DY oc A 2 . In this letter, we 
use a measured b DY sa 0.022 GeV 2 to fix A 2 and 
j^sh Taking the small xdy limit, and using the fact 
that Rj) Y (Q 2 ) w 1, we derive 



otDYiqr) ~ 1 



>DY 
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In deriving Eq. (jioj) , we used A 1 ! 3 ~ hi(A), which is a 
good approximation for most relevant targets. Eq. ( p^| ) 
shows that the leading contribution to aDY(qr) does not 
depend on the A. 
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where (q 2 ) hA = (g 2 ) hN +A(q 2 , ) hA . A{q 2 T ) hA is the trans- 
verse momentum broadening and calculable in QCD per- 
turbation theory In Eq. (|), T{q%) hA is a small 
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FIG. 1. dDy(9r) for the Drell-Yan production as a func- 
tion of qr- At small qr, the aoYiqr) in Eq. ^ is used. At 
large qr, QCD predictions from Ref. JTJJ] are plotted with 
ro = 1.1 fm and [xG(x)] xx o ~ 3 for the possible maximum 
and minimum values of the C. 
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In Fig. [j], we plot ccdy(<1t) for targets with A = 12 
and A = 184. At small qx, we used Eq. (||) and 
set the small T(q^) — 0. We also used the measured 
= 1 - 8 Gey2 ' and A (<lT} hA = 0.022A 1 /3 QeV 2 
Jl2|,[r7| . At large qx, we plot the QCD predictions from 
Ref J10|, which depend on both T SH and T DH . T SH is 
fixed by A(q T ) hA . However, there is no direct observable 
yet to extract T DH [jllj. Because of the operator defini- 
tion of T DH (x%, X2), it was assumed that T^ 2 (x%,X2) — 

(2ttC)A 4 / 3 fi(xi) fi{x2) with fx and f% parton distribu- 
tions of flavor fi and / 2 ficfl - Assuming no quantum 
interference between different nucleon states, one derives 
C = 0.35/(87rr 2 ) GeV 2 with r w 1.1 - 1.25 fm, which 
is just a geometric factor for finding two nucleons at the 
same impact parameter jl(J. On the other hand, when 
x\ (or X2) goes to zero, the corresponding parton fields 
reach the saturation region, and the T DH is reduced to 
T SH . Therefore, we have C « A 2 / (2n[xG(xY\ x ~ ) ©, 
where [xG(x)] x ~q is of the order of unity |Q9[. Because 
of a combination of a small value of the measured A 2 
from Drell-Yan data |l7]] and a choice of [xG(x)] x ~o « 3 
fL2"f , these two approaches result into a factor of 20 dif- 
ference in numerical value for the parameter C ]l2] | . The 
value for the C obtained in Ref. [ |10f without any quan- 
tum interference should represent a possible maximum 
for C, while the value obtained in Ref. jl2| with full 
quantum interference (in saturation region) represents 
a possible minimum for the C . In Fig. [l], we plot the 
perturbatively calculated aDY(qr) fl0| at large qx with 
the maximum and minimum values of the C discussed 
above, and let the a(qx) in small qx naturally linked 
to that at large qx. Fig. [I] also shows that aDY(qr) 
is insensitive to the atomic number A. In Fig. we 
plot R(A,q T ) = A aDY ^ by using aDY(qr) in Figjl] 
and compare our predictions with data from E772 |p|,]l8P 
Without any extra free fitting parameters, our predic- 
tions shown in Fig. || are consistent with data at small 
qx, and due to the large error in data at high qx, current 
Drell-Yan data are consistent with almost any value for 
the C between the maximum and the 
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FIG. 2. Comparison of our cioy((1t) in Eq. (g|) with the 
Drell-Yan data from Fermilab E772 collaboration [|,^J) . 

minimum discussed above. In order to test the theory 
and pin down the value of the C , we need either better 
data or different observables. Because of the strong de- 
pendence on the double-hard subprocesses, the angular 
dependence of the Drell-Yan pair jll]] could be an excel- 
lent signal for measuring the C. 

Similarly, we can also obtain the a{qx) for J/ip pro- 
duction. Kinematically, hadronic J/ip production is 
like Drell-Yan production with Q ~ Mj , and its qx- 
spectrum can also be characterized by three different re- 
gions as in the Drell-Yan case. Because J/ip mass Mj/^ 
is smaller than any typical Q measured for the Drell- 
Yan continuum, the logarithm (a s ln 2 (Mj/^/q T ) ) n for 
J/ijj production is less important. Consequently, at fixed 
target energies, a Gaussian-like distribution can fit J/^'s 
<7r-spectrum even better. Therefore, following the same 
arguments used above for the Drell-Yan production, we 
derive a(qx) for 3/ip production in small qx region, 
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where Rfo = (l/A)afy/a™ 



'/</> 



and Xj/?p is defined by 
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Similar to the Drell-Yan case, r(g^)W and AF(q^)W^ 
are perturbatively calculable, and much smaller than 
(q T ) h 3 ^ and A(<z 2 )^, respectively @. 

One major difference between J/^p production and the 
Drell-Yan process is the nuclear dependence of Rf/^- 
Clear nuclear suppression for (Tj^ has been observed Q . 
Since we are only interested in the general features of 
a(qx), we adopt the following simple parameterization 
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which fits all experimental data on 3/tp suppression in 
hadron- nucleus collisions |2]]J23]. The A 1 / 3 factor in 
Eq. (|l^) represents an effective medium length, and the 
(3 is a constant [[12|. 



Similar to the Drell-Yan process, A(q T )'j A ^ is propor- 
tional to four-parton correlation functions [ |12| . Due to 
final-state interactions for 3/ijj production, A(q T ) l ^ A ^ de- 
pend on both quark-gluon and gluon-gluon correlation 

b m AV\ 



functions. We can define [|12| that A(q T ) hA 



where bj/^p can be calculated or extracted from data. 
From Ref. 0, we obtain b 3Jil w 0.06 GeV 2 . With 
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(<?t)jM ~ L68 GeV2 0! we Can als0 take the Sma11 

X] m limit in Eq. (|i~l| ) and obtain 



12 \h 

Wt/j 
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where A 1 / 3 ~ ln(A) was again used. It is clear from 
Eq. ( fl4"| ) that in small qx region, aj/^(qx) is also in- 
sensitive to the atomic number A of targets. Further- 
more, the nuclear suppression in Rf^ corresponds to a 
(^-independent shift in the magnitude of aj/^,(gr). 
Because the Rf/^ , Xj/V> i and other physical quantities 



in Eq. ( |ll| ) can depend on if, otj/^,(qx) can also be a 
function of xp. Experiments show that the larger xf, the 
more suppression for J ftp production (or smaller Rf^) 



J23[ . Consequently, from Eq. <]1 1|), we will have smaller 
aj/^(qx) at larger xp, which is consistent with exper- 
imental data ||. Although we do not have all needed 
physical quantities in Eq. (|Il]) for predicting the aj in 
different xf regions, we can still test the universality of 
ctj/^(qx)'- the quadratic dependence on qx- In Fig. ||, we 
plot our fits using OLj/^qx) in Eq. (^lj) and compared it 
with E866 data in three xp regions: small (SXF), inter- 
mediate (IXF), and large (LXF). Clearly, our universal 
functional form for otji^qx) is consistent with all data 
in small qx region (qx < q T ~ Mj/^,/2). 

In summary, we derived an universal functional form 
for a(qx) for both the Drell-Yan and production in 
small qx region (qx < q^ = kQ with k ~ 1/3 — 1/2). 
All parameters defining a(qx) in Eqs. (||) and ( p"l| ) are 
completely determined by either perturbatively calcula- 
ble or independently measurable quantities. We show 
that a(qx) is extremely insensitive to the atomic weight 
A of targets. For the Drell-Yan process, a(qx) in Eq. (||) 
can be naturally connected to the perturbatively calcu- 
lated a{qx) at large qx pX( |. A similar test can also be 
carried out for J/?p production. suppression in rel- 
ativistic heavy ion collisions was predicted to signal the 
color deconfinement p4| . On the other hand, significant 
3/tp suppression has been observed in hadron-nucleus col- 
lisions pQ] . Therefore, understanding the features ob- 
served in a(qx) for J/ip production is very valuable for 
finding the true mechanism of J/ip suppression. 
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FIG. 3. Comparison of our fits, using the Qj/^,((?t) in 
Eq. ([n]), with J/i/j data from Fermilab E866 collaboration 
H in three different regions of xf- 
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